Bound states of dipolar molecules studied with the Berggren expansion method 
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Bound states of dipole-bound negative anions are studied by using a non-adiabatic pseudopotential 
method and the Berggren expansion involving bound states, decaying resonant states, and non- 
resonant scattering continuum. The method is benchmarked by using the traditional technique 
of direct integration of coupled channel equations. A good agreement between the two methods 
has been found for well-bound states. For weakly-bound subthreshold states with binding energies 
comparable with rotational energies of the anion, the direct integration approach breaks down and 
the Berggren expansion method becomes the tool of choice. 
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I. INTRODUCTION 

Weakly-bound many-body systems are intensely stud- 
ied in dilferent domains of mesoscopic physics [HQ, in- 
cluding nuclear 0-01' molecular (sl-lisj. and atomic (TJ- 
[TgI physics. In this context, dipolar anions are one of the 
most spectacular examples of marginally bound quantum 
systems (TMll- 

The mechanism for forming anion states by the long- 
range dipolar potential has been proposed by Fermi and 
Teller [301 , who studied the capture of negatively charged 
mesons in matter. They found that if a negative meson 
is captured by a hydrogen nucleus, the binding energy of 
the electron becomes zero for the electric dipole moment 
of a meson-proton system ficr = 1-625 D. Later this result 
was generalized to the case of an extended dipole with an 
infinite moment of inertia [40^ . Lifting the adiabatic ap- 
proximation by considering the rotational degrees of free- 
dom of the anion (T7l-[22] turned out to be crucial; it also 
boosted the critical value of /i to about 2.5 D. For anions 
with fj, > ficr, the number of bound states of the electron 
becomes finite, and the critical electric dipole moment 
/icr depends on the moment of inertia of the molecule. 
In the non-adiabatic calculations, the pseudo-potential 
was used to take into account finite size effects, repulsive 
core, polarization effects, and quadrupolar interaction. 
The pseudo-potential method has provided a convenient 
description of binding energy of the electron bound by an 
electric dipolar field. Recently, this method was applied 
to linear electric quadrupole systems 



41|. Some recent 



theoretical studies of dipole-bound anions also employed 
the coupled cluster technique [42-44]. 

The unbound part of the spectrum of multipolar an- 
ions has been discussed theoretically in Refs. [i^ HI] 
and Refs. quoted therein. Resonance energies of 
dipolar anions have been determined experimentally by 
low energy electron scattering off the dipolar molecules 



[2511261,1291,1311,1331. 

Both the long-range dipole potential and the weak 
binding of dipolar anions provide a considerable chal- 
lenge for theory. The impact of the molecular rotation on 
a weakly-bound electron can be represented by coupled- 
channel (CC) equations that can be solved by means of 
the direct integration. While this approach correctly pre- 
dicts the number of bound states of polar anions, it is 
less precise for treatment of weakly-bound excited states. 
Moreover, it cannot be used for studies of dipolar anion 
resonances because the exact asymptotics for a dipolar 
potential in the presence of a molecular rotor cannot be 
determined. 

In this paper, we apply the complex-energy configura- 
tion interaction framework based on the Berggren ensem- 
ble [13] to the problem of bound states in dipole-bound 
negative anions. The Berggren completeness relation is a 
resonant-state expansion; it treats the resonant and scat- 
tering states on the same footing as bound states. We 
have successfully applied this tool to a variety of nuclear 
structure problems pertaining to weakly-bound and un- 
bound nuclear states [isl - lsij (for a recent review see Ref. 
(5^). The nuclear many-body realisation of the complex 
energy configuration interaction method is known under 
the name of the Gamow Shell Model. 

Resonances do not belong to the Hilbert space, so 
the mathematical apparatus of quantum mechanics in 
Hilbert space is inadequate for Gamow states [s^ , which 
are not square-integrable. It turned out that the mathe- 
matical structure of the Rigged Hilbert Space (RHS) [s^ - 
[56| can accommodate time-asymmetric processes, such as 
particle decays, by extending the domain of quantum me- 
chanics. The mathematical setting of the resonant state 
expansions follows directly from the formulation of quan- 
tum mechanics in the RHS [13, , rather than the usual 
Hilbert space [H-IH. 

The Berggren ensemble provides a natural generaliza- 
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tion of the configuration interaction for the description 
of the particle continuum. The complex-energy Gamow- 
Siegert states [H, [s^ states have been used in various 
contexts in nuclear, atomic, and molecular physics (60l - 
[73} . Some recent applications of Gamow-Siegert states, 
also in the context of a CC formalism relevant to the 
problem of dipole anions, can be found in, e.g., Refs. 

This paper is organized as follows. The Hamiltonian 
of the pseudo-potential method is briefly discussed in 
Sec. ini The CC formulation of the Schrodinger equation 
for dipole-bound negative anions is outlined in Sec. IIIII 
Section lTVl discusses the direct integration method (DIM) 
for solving the CC problem with a focus on difficul- 
ties in imposing proper boundary conditions when the 
rotational motion of the molecule is considered. The 
Bergggren expansion method (BEM) is introduced in 
Sec. |Vl Section IVII specifies the coupling constants of 
the pseudo-potential and other calculation parameters. 
Salient features of DIM and BEM solutions are com- 
pared in Sec. IVlIl The predictions of DIM and BEM for 
low-lying energy states and r.m.s. radii of Lil~, LiCl~, 
LiF~, and LiH~ anions are collected in Sec. I Villi Finally, 
Sec. IIXI contains the conclusions and outlook. 



rUe its mass. The interaction V is approximated by a 
one-body pseudo-potential V{r, 6) acting on the valence 
electron [^[7i,|80|: 

V{r,e) = V^{r,e) + V^{r,e) + VQ^Sr,0) + VsK{r), (3) 

where 9 is the angle between the dipolar charge separa- 
tion s and electron coordinate; 

V,{r,e) = -^ie f?)';^^A(^o«^) (4) 

A=l,3, 



is the dipole potential of the molecule; 



K.(r,0) --^[ao + a2P2(cos^)]/(r) 



(5) 



is the induced dipole potential, where ao and 02 are 
the spherical and quadrupole polarizabilities of the linear 
molecule; 



-^QzzP2icos9)f{r) 



(6) 



is the potential due to the permanent quadrupole mo- 
ment of the molecule; and a short-range potential 



II. HAMILTONIAN 



VsR{r) = Vbexp(-(r/rc) ) 



(7) 



A dipole-bound negative anion is composed of a neutral 
polar molecule with a dipole moment greater than ficr 
and a valence electron. The Hamiltonian of the total 
system can be written as: 



Htot — He+ H,nol + V 



(1) 



where He is the Hamiltonian of the valence elec- 
tron, Hmoi is the Hamiltonian of the molecule, and V 
is the electron-molecule interaction. The many-body 
Schrodinger equation for Htot couples all electrons of the 
system; hence, an approximation scheme has to be de- 
veloped. 

As a first simplification, we assume that the vibrational 
motion of a molecule is much slower than other modes so 
that it can be treated in the Born-Oppenheimer approx- 
imation. The Hamiltonian ([1]) simplifies considerably if 
one considers anions of closed-shell systems. Moreover, 
if spin is neglected 22] , the molecule can be treated as a 
rigid rotor. Note that the energy scales associated with 
the rotational motion of the molecule and the motion of 
the weakly-bound valence electron may be comparable. 
Consequently, there appears a strong non-adiabatic cou- 
pling between the molecular angular momentum j and 
the orbital angular momentum £ of the electron. Eq. ([T]) 
thus writes within this approximation scheme: 
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(2) 



where / is the moment of inertia of the neutral molecule, 
Pe is the linear momentum of the valence electron and 



accounts for the exchange effects and compensates for 
spurious effects induced by the cut-off function 



/(r) = 1 - exp{-(r/ro) } 



(8) 



introduced in Eqs. (I5I6P to avoid a singularity at r — 0. 
The parameter rg in Eq. ([5]) is an effective short-range 
cutoff distance for the long-range interactions. 



III. COUPLED-CHANNEL EXPRESSION OF 
THE HAMILTONIAN 

The eigenfunctions of the Hamiltonian ([2]) can be con- 
veniently expressed in the CC representation: 



(9) 



where the index c labels the channel, u'l{r) is the radial 
wave function of the valence electron in a channel c, and 
the channel function $^ ^ arises from the coupling of jc 
and to the total angular momentum J of the anion: 
j + £ = J . Due to rotational invariance of Htot, its ma- 
trix elements arc independent of the magnetic quantum 
number M , which will be omitted in the following. 

The potential V{r,6) in Eqs. (I3]-[7|) can be expanded 
in multipoles: 



y(r,0) =^yA(r)PA(cos0), 



(10) 
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where 



(11) 



The matrix eleraents of P\{cos9) between the channels 
c and c' are obtained by means of the standard angular 
momentum algebra: 



(<i>/^,,, I Pa (cos 0)1$/^, J 



jc A jc 




Ic' A 





X V(24' + 1)(24 + l)(2jc' + l){2jc + 1). 



(12) 



In the following, we express r in units of the Bohr ra- 
dius oo, / in units of meOg, and energy in Ry. The radial 
functions (?") are solutions of the set of CC equations: 



4(4 + 1) ]c{3c + 1) 



E, 



<{r) 



(13) 



where Ej is the energy of the system and 



viAr) = V ($/ „ IPa(cos^)|$'^, ) Vy.{r). (14) 



A 



IV. DIRECT INTEGRATION OF 
COUPLED-CHANNEL EQUATIONS 



The CC equations (|T3l) can be solved by the DI 
method. Below we describe the method used to generate 
the channel wave functions Uc{r) (from now on, the quan- 
tum number J is omitted to simplify notation) obeying 
the physical boundary conditions. Namely, we assume 
that Uc{r) is regular at origin: Uc{r = 0) = 0, and for 
r — )• +00 it behaves like an outgoing wave u^{r). 

The central issue of DI lies in the boundary condition 
at infinity. Indeed, as we shall see in Sec. (jIV Bp . an 
asymptotic wave function of a dipole-bound anion is not 
analytic in general, so that one cannot exactly impose 
outgoing boundary conditions. This calls for the use of 
controlled approximations. In the following, we describe 
the numerical integration of CC equations. While the 
method is standard (cf. Sec. 3.3.2 of Ref. [HI), this 
particular application is not; hence key details should be 
given. 



A. The basis method with the direct integration 

To integrate CC equations, we introduce the matching 
radius that defines the internal region [0 : r™], where 
the centrifugal potential is appreciable, and the external 
zone [fm : -|-oo]. An internal basis function u^.^(r) in 



[0 : r™] is regular at r = 0: 



<c,(0 



(15) 



in one channel Cf, = {jcbAc^)- The CC equations imply 
that when r — >■ the internal channel wave functions 



''b-c 



ir) with c ^ cjj must behave as: 



2meVccM 



ln(r/r„i) 



(«c,+2)(fe,+3)~^e(<'c + l) 



for 4 = 4, + 2, 
otherwise. 



Note that is it necessary to pay attention when integrat- 
ing CC equations close to r — s, as the potential ^ is 
not differentiable therein. 

In the external region [r„i : -l-oo], the basis wave func- 
tions are denoted u^.f.^ {r) . By construction, at very large 
distances of the order of hundreds of (asymptotic re- 
gion), u+ (r) ^ for Cb = (jc,,4J and ■ti^^(r) = for 



other channels c ^ Cb- The asymptotic behavior of exter- 
nal channel functions is discussed in Sec. IIV Bl below. 

Both sets of internal and external basis functions are 
used to expand the channel fimction Uc{r): 



Uc{r) 



bjr) 



for r < r„ 
for r > r.„ 



The matching conditions at r — 

J2 [C'b<cirm) - C+u+Arra) 



E 



CI 



dr 



dul 

-(^m) — Cfc {fm) 



= 0, 



0, 



(16) 

(17) 
(18) 



form a linear system of equations: AX = 0. The condi- 
tion of det A = determines the energy of a bound or res- 
onant state. (One can thus see that det A is thus the gen- 
eralization of the Jost function for CC equations.) Once 
the eigenenergy has been found, the amplitudes C°,C^ 
are given by the eigenvector X oi A. The overall norm 
is determined by the condition: 



J2 K{r)\'dr^l. 



(19) 



B. The coupled-channel equations in the 
asymptotic region 

At large distances, Vcc' (r) can be written as: 



Vcc'ir) 



2ni„ 



^3(0 



(20) 



where Xcc' is a constant and V^{r) decreases for r — !■ +oo 
as . In the following, we shall assume that V^{r) = 
in the asymptotic region. As the numerical integration up 
to r ~ 100 ao is stable, the error made by neglecting V3 
is around 10~^ olq^ ^ which is sufficiently small to insure 
that the asymptotic zone has been practically reached. 

Let us first consider the case of an infinite moment of 
inertia / — ^ +00. Here, Eq. (fT3|) becomes: 



<{r) 



4(4 + 1) 



{r)+Y,^^Ar)-k^u,{r), (21) 
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where k = s/E. The outgoing solution of (|2ip in a basis 
channel h can be written in terms of spherical Hankel 
functions: 



(22) 



where ^^j^ is an effective angular momentum given by 
eigenvalues of the eigenproblem 



Indeed, it immediately follows from Eqs. (PT|) and 
that: 

<c(0" = { "^MSllA _ ,A ,+^(,), (24) 



complete basis of single-particle states. Since our goal 
is to describe weakly-bound or unbound states, special 
care should be taken to treat the asymptotic part of 
wave functions as precisely as possible. A suitable ba- 
sis for this problem is the one-body Berggren ensemble 
[47I l63l . Is^ . This basis is generated by a finite-depth 
spherical potential and contains bound (6), decaying (c?), 
and scattering (s) one-body states. Fot that reason, the 
Berggren ensemble is ideally suited to deal with struc- 
tures having large spatial extensions (such as halos or 
Rydberg states) or outgoing behavior (such as decaying 
resonances). Some recent applications, in a many-body 
context, have been reviewed in Ref. l52t. 



A. The Berggren basis 



in terms of an ef- 



so the physical interpretation of v^^ ^ 
fective angular momentum is justified. 

If / is finite, however, solutions of Eq. are no 

longer analytical at large distances. Nevertheless, it 
is possible to construct an adiabatic approximation for 
Uc(r) in the asymptotic region. To this end, one de- 
fines the linear momentum fcc = \JE — jc{jc + 1)/^ for a 
channel c. In the asymptotic region, Eq. (jl3l) becomes: 



<'(r) = 



^u,{r)+y^^Uc'ir)-klucir), (25) 



where, compared to Eq. (PT|) . k is replaced by the chan- 
nel momentum kc- This approximation can be applied 
if \E\ ^ jc{jc + l)/-f for all channels of importance. In 
those cases, one can introduce an ansatz for Uc{r) by re- 
placing k by kc in Eq. (|22p . The relative error on a basis 
function M^''^''(r) associated with this approximation is 



E 



Xcc 9, 



'■iff 



- 1 



(26) 



i.e., is of the order of |fcc — fcc'|/r^. 

In practical calculations, / ^ 10''' and jmax ~ 7. 
This gives jmax{jmax ^ / 1 ^ 10^^ ■ Consequently, if 
\E\ > 10-3 Ry, the error % - kc'\/r^ < 10~^ Oq^ for 
r ~ 100 oq is close to that associated with the neglect of 
V3 (r) . On can thus see that the proposed ansatz accounts 
for the coupling term (|20p in many cases. However, this 
approximation breaks down for weakly-bound/unbound 
states with \E\ < 10"^ Ry; hence, a more adequate theo- 
retical method based on a resonant state expansion needs 
to be introduced. 



V. DIAGONALIZATION WITH THE 
BERGGREN BASIS 

Another way to find eigenstates of the CC problem 
is to diagonalize the associated Hamiltonian in a 



The finite-depth potential generating the Berggren en- 
semble can be chosen arbitrarily. To improve the conver- 
gence, however, it is convenient in practical applications 
to use a one-body potential, which is as close as possible 
to the Hartree-Fock field of the Hamiltonian in question. 
Therefore, in the case of the one-body problem (|13p . the 
most optimal potential to generate the Berggren basis is 
the diagonal part of Vcc' (r) ■ This means that the basis 
states ^k.cir) are eigenstates of the spherical potential 
Vccir): 

Kcir) = (^^^^^ + vUr) - fc^) ^kAr) (27) 
that obey the following boundary conditions: 



$fc,c(r) - C"r 



for all states. 



(28) 

$fc,c(r) - C+H+ikr), r ^ +00 (6,d), (29) 

$fc,c(r) - C+H+{kr) + C-Hl^{kr), r ^ +00 (s), (30) 

where the boundary conditions at r ~ and at r — !■ +00 
for scattering states (s) are standard, and for bound and 
decaying states (6, d) one imposes the outgoing boundary 
condition. Note that k in Eq. (P?]) is in general complex. 

The scattering states are normalized to the Dirac delta, 
which results in a condition for the C~ and ampli- 
tudes in ([50)1 [ivt: 



1. 



(31) 



The normalization of bound states is standard as well, 
but that of decaying resonant states is not. Indeed, reso- 
nant states rapidly oscillate and diverge exponentially in 
modulus along the real r-axis; hence, one cannot calcu- 
late their norm in the same way as for the bound states. 

The solution of this problem is provided by the exterior 
complex scaling |83(|, i.e., one calculates the norm of the 
resonant state using complex r radii: 



+ 00 



$fe^(i? + a;e'^)e*''da;,, (32) 
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where i? is a radius taken sufEciently large so that con- 
dition (|29l) is fulfiUed. In the above formula, is an 
angle of rotation chosen so that ^k,c{R + xe^^) for 
X — >■ +00, which is always possible provided 9 is larger 
than a critical value depending on k [s^l- Note that no 
modulus enters Eq. ([32l) . This arises from the finite- 
lifetime character of resonant states, which requires us to 
use the biorthogonal scalar product J2, 84- 87] . It can be 
shown [sB - lsrj that the norm defined in Eq. (l35t is indeed 
independent of R and 9, as expected from a norm. Since 
the expression ([5^ is also valid for bound states, bound 
and decaying states enter the Berggren ensemble as one 
family of resonant states. 

The exterior complex scaling can be used to calculate 
matrix elements of a one-body operator 0{r) as well, 
provided it decreases faster than 1/r along the complex 
r-contour: 

{^k'.c'\0\^k,c)^ f ^k',c'{r) Occ'ir) <^k,c{r)dr 

JQ 

r+oo 

+ ^k'.c'{zix))Occ'{zix))^k,M^))e''dx, (33) 
Jo 

where z{x) — R + xe^^, and \^k,c) and \^k',c'} can here 
be bound, decaying, or scattering states. 

B. The Berggren completeness relation 
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FIG. 1. Berggren ensemble in the complex momentum plane. 
The bound (b) and antibound (a) states are distributed along 
the imaginary fc-axis at Im{k) > and Im{k) < 0, respec- 
tively. The decaying resonant states (d) are located in the 
fourth quadrant {Re{k) > 0, /m(fc) < 0. The Berggren com- 
pleteness relation involves bound states, scattering states (s) 
on the contour, and decaying states lying between the 
real-fe axis and If antibound states are included, the £^ 
contour has to be slightly deformed [s^. [s^. If the contour 
lies on the real fc-axis, the Berggren completeness relation 
reduces to the Newton completeness relation ^,,8^ involving 
bound and real-energy scattering states. 

Figure [T] shows a distribution the Berggren ensemble 
in the complex momentum plane. To determine the ba- 
sis, one first chooses a £^ contour in the fourth quad- 



rant containing the decaying eigenstates. The scattering 
states of the ensemble lie on this contour. The resonant 
part of the ensemble contains the bound states lying on 
the imaginary- fc axis and those decaying states of ^7} 
that are found in the region between the real-fc axis and 
£+. The Berggren basis is built from all those states: 

J2 l*fcn,c>(*fc„,c|+ / |*fc,c)($fe,c| dk^l. (34) 

ne{b,d) 

This completeness relation corresponds to a given chan- 
nel c; hence, one has to construct Berggren ensembles for 
all the channels considered in Eq. (|13p . 

In order to be able to use (IMl) in practice, one needs 
to discretize £+. Our method of choice is to apply 
the Gauss-Legendre quadrature to each of the segments 
defining £+ in Fig. [T] The last segment, chosen along 
the real-fc axis, extends to the large cutoff momentum 
fc = kmax that is sufficiently large to guarantee complete- 
ness to desired precision. It is then convenient to renor- 
malize scattering states using the corresponding Gauss- 
Legendre weights ujkn ■ 

|*n,c) = ^A:^|$fc„,c). (35) 

The discretized Berggren completeness relation, used in 
practical computations, reads: 

N 

5^1$,,,) ($,,,1-1, (36) 

4=1 

where the N basis states |<I>,;.c) include all bound, decay- 
ing, and discretized scattering states of the channel c. By 
using Eq. ([55l) . the Dirac delta normalization of scatter- 
ing states has been replaced by the usual normalization 
to Kronecker's delta, in this way, all \^i.c) states can be 
treated on the same footing in Eq. p6p . as in any basis 
of discrete states. 



C. Hamiltonian matrix in the Berggren basis 

As the basis states |$i,c) are generated by Vcdr), the 
Hamiltonian matrix within the same channel c is diago- 
nal: 

($.',c|/»|<I>.,c) = (fcf + MZl±l)^ Sw (37) 

Matrix elements between two basis states belonging to 
different channels c and c' are: 

Jo 

+ / <^^',c'{z{x)) Vcc'{z{x))^^Az{x))e'^ dx, (38) 
Jo 
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where the complex scaling ((33t can be used, because 
Vcc'{r) decreases at least as fast as r~^. 

As the off-diagonal matrix elements are present only 
for c 7^ c', the Berggren basis generated by Eq. (f27|) is 
optimal. The channel wave functions Uc{r) can be ex- 
pressed in the Berggren basis by diagonalizing the matrix 
of /i (I37l38)l . 



VI. CALCULATION PARAMETERS 

Results of the direct integration method (DIM) depend 
both on the parameters of the pseudo-potential ^ and 
on the cutoff value of the electron orbital angular mo- 
mentum imax considered in the CC problem. They are 
fixed to reproduce the experimental value of the ground 



-4.483- 10"^ Ry 



state energy of the LiCl an 
27j. 

The most important term in ([3|) is the dipole poten- 
tial Vfj,, which depends only on the dipole moment fi 
and the size s of the neutral molecule. The remain- 
ing parameters of the pseudo-potential are taken from 
Ref. [2^, namely: = 15.3 aj], 012 = 1.1 09, fo = 2.2 oq, 
rc = 2.828 flo, Q,, = 3.28 eag, and Vq = 2.0 Ry. The 
moment of inertia parameters are: / = 150, OOOmeflo 
for LiCr, 240, 000 meag for Lil", 82,000meag for LiF", 
and 26,000TOeag for LiH~. The dipole moment of each 
molecule considered in this work is known experimentally 
and has been taken from the NIST database. 
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FIG. 2. 

near r = for the first excited = Oj state of LiCP 
calculated in DIM (solid line) and BEM (dotted line) with 
imax = 9. The charge separation s of LiCl has been adjusted 
in both approaches to the experimental ground state energy 
in the limit (.max — >■ 00. 

For imax = 9 the ground state energy of the LiCl^ 
anion is reproduced by taking the charge separation 
0.336 cq. To remove the dependence of results 



(9) 
■^DIM 

on £rn 



extrapolated for imax — > 00, and the size of the charge 
separation s is adjusted to reproduce the experimental 
binding energy. In this case, s[^j^ — 0.337 oq. The 
matching radius was taken as — a-o- This value was 
found to optimize the DIM procedure. 

Anion spectra in the BEM depend sensitively on the 
cutoff parameter kmax of the single-particle basis. How- 
ever, as we shall see in Sec. lVIIl for a chosen value of kmax 
they are practically independent of imax- In this study. 



we have chosen k„ 



1.53 Oq for each partial wave in 



in the DIM, the ground state energy of LiCl is 



order to attain both a good numerical precision and ap- 
proximately the same value of the dipole size parameter s 
as in DI. In this case, Sggjy; — s'^l^ = 0.336 oq. We have 
used complex contours with straight segments connecting 
points: fci = (0,0), /ca = (0.15, -i0.04), ka = (1,0), and 
^4 = kmax m units of Oq ^. Each scattering contour has 
been discretized with 220 points. The precise form of the 
contour does not change results; since the applications 
carried out in this work pertain to bound states only, we 
could have used real scattering contours, i.e., the Newton 
completeness relation [H, . 



VII. NUMERICAL TESTS AND 
BENCHMARKING 

Along with the asymptotic behavior of channel wave 
functions, treated approximately with the DIM and ex- 
actly within BEM, the Hamiltonian 1^ cannot be iden- 
tically represented in both approaches. Indeed, since 
the potential V'^(r) ^ is not differcntiable at r = s, it 
cannot be treated exactly in BEM because the channel 
wave functions expanded in the Berggren basis are ana- 
lytic by construction. In practice, this translates into a 
node beyond r = s in DIM channel wave functions, which 
is absent in BEM. This is illustrated in Fig. [2] for a 
(j = 0, ^ = 0) channel function corresponding to the first 
excited = 0^ state of LiCl~. It is to be noted, how- 
ever, that beyond this point the channel wave functions 
calculated with both methods are very close and - as will 
be discussed later - this near-origin pathology has a very 
small impact on the total energy as the contribution from 
this region is small. 

As discussed in Sec. lIVBi DIM is inadequate for states 
with very small energies, while BEM has been shown 
to be very precise in this case. On the other hand, for 
states with binding energies typically greater than 10~^ 
Ry, BEM yields channel wave functions that exhibit spu- 
rious low-amplitude oscillations. Figure [3] illustrates such 
wiggles in the tail of the channel wave function Mj=o,f=o 
of the — 0^ ground state of LiCl^. For such well- 
bound states that quickly decay with r, the standard size 
of the Berggren basis (measured in terms of contour dis- 
cretization points and kmax) is not sufficient. The DIM 
is thus preferable for such cases, as the asymptotic be- 
havior of well-bound states is treated almost exactly (see 
Sec. HVBl) . 

The direct integration becomes numerically unstable 
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100 
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FIG. 3. The modulus of the channel wave function Uj=o,<=o 
for the = 0^ ground state of LiCl~ calculated in DIM 
(solid line) and BEM (dotted line) with Imax = 9- At large 
distances, spurious wiggles appear in BEM results (see the 
inset) due to basis truncation. 



when the channel orbital angular momentum becomes 
large, around £c — 10, even for the states with relatively 
large binding energies. In this case, the matrix of basis 
channel wave functions u'^-ci^m) and u'^.^{rm) and their 
derivatives, introduced in Sec. IIV Al in the context of 



matching conditions at r 



is ill-conditioned and its 



eigenvector of zero eigenvalue becomes imprecise. This 
results in a discontinuity at and spurious occupa- 
tion of channels with large orbital angular momentum 
4 > 10. This is iUustrated in Fig. H for the ,r = Oj^ 
ground state of LiCl^. A a result, the energy and spa- 
tial extension of the electron cloud distribution of the CC 
eigenstate become incorrect. 

The convergence of the LiCl^ ground state energy with 
respect to £max is shown in Fig. [5] One may notice 
an exponential convergence of calculated DIM energies 
with £max for 6 < imax < 10 and a clear deviation for 
£ > 11, which is related to the discontinuity of channel 
wave functions for £c > 10. The energy calculated in 
BEM is perfectly stable with imax- 

The rapid converge of BEM with i is due to kmax- 
truncation of the single-particle basis that suppresses 
contributions from large-£ configurations. This is illus- 
trated in Fig. [51 which displays the average modulus of 
the off-diagonal matrix element of the channel-channel 
coupling in BEM: 



1 ^ 



(39) 



n.n' 



between the first channel c — {j = 0,£ ^ 0) and higher- 
£ channels c'. Only the channels with < 5 and 
\£c — £c'\ < 3 contribute significantly to the channel cou- 
pling matrix element. We checked that this is generally 
the case. Using the same truncation, the DIM yields 



o 
:3 




0.2 0.4 0.6 0.8 1.0 1.2 

r(ao) 

FIG. 4. The modulus of the channel wave function Uj^o,e=o 
for the = 0^ ground state of LiCl~ calculated in DIM with 
several values of imax- For imax > 10, one may notice the 
development of a discontinuity at the matching point r,„ = 
ao. In such cases, the channel wave function becomes ill- 
conditioned, introducing serious errors in CC eigenenergy and 
eigenfunction. 



-0.044543 



0.044677 



-0.044812 




7 8 9 10 11 12 

^max 



FIG. 5. The dependence of the LiCl ground state energy on 
imax in DIM (dots) and BEM (triangles). The DIM results 
converge exponentially (red line) . This allows us to determine 
the asymptotic value of energy at £max oo. 



numerically stable results. In this case, the energies of 
well-bound states {\E\ > 10"^ j^y) 

agree in both meth- 
ods. 

The numerical instability of DIM at large £max leads 
to a collapse of calculated radii. Figure [7] shows the de- 
pendence of the ground state r.m.s. radius of LiCl~ on 
£max- This result, together with discussion of Fig. [51 sug- 
gests that the BEM can provide practical guidance on the 
minimal number of channels in the CC approach. 

In practical applications, spurious oscillations in BEM 
channel wave functions for well-bound states can be taken 
care of by extrapolating wave functions from the inter- 
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FIG. 6. Average off-diagonal matrix element Ao,c (|39p of the 
channel-channel coupling in BEM between the channel (j = 
0,^ = 0) and c' for the J"" = 0+ ground state of LiCl". 



BEM 



LiCI- 

y"=o; 



DIM 



10 



11 12 



FIG. 7. The dependence of the LiCl ground state r.m.s. 
radius on £^ax DIM (dots) and BEM (dotted line). The DIM 
results are stable up to imax = 10. 



mediate region of r, where they are reliably calculated, 
into the asymptotic region. This can be done by applying 
the analytical expression: 



Uc{r) 



lini Ur(r) 



M 



J 



(40) 



(c) 

where kc is the channel momentum and a - are param- 
eters to be determined by the fit. The precision of this 
procedure can be assessed by computing tlie norm of the 
eigenstate. Using this procedure, one obtains perfectly 
stable r.m.s. radii in BEM for different values of tmax, 
as can be seen in Fig. [T] 

Figures I51I10I compare the four most important chan- 
nel wave functions (£, j) of DIM and BEM corresponding 
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FIG. 8. Most important channel wave functions Uc(r) with 
c = ( j, t) for the — 0^ ground state of LiCl~ , as calculated 
in DIM (solid hue) and BEM (dashed line) with tmax = 9. 




FIG. 9. Similar to Fig. [8] but for the first excited J" = Oj 
state of LiCl^ . 



to the three lowest — 0^ eigenstates of LiCl~. For 
the ground state, both approaches predict the same en- 
ergy E = —4.483 • 10"^ Ry and the channel functions 
are practically identical. For the first excited state, the 
agreement is still reasonable. Here, the energy in DIM 
is = -7.374 • 10-* Ry while BEM gives slightly more 
binding: E = -8.241 • lO^^Ry. Consequently, the BEM 
wave functions decay faster than those computed with 
DIM. For a second excited Qt, state, both methods differ 
markedly. This state has a sub-threshold nature, with 
-Bdim = -7.051 • 10-6 Ry and Sbem = -9.907 • 10"^ Ry. 
For this extremely diffused state, the direct integration 
method fails completely. This is manifested by the very 
different nodal structure of channel wave functions in 
DIM seen in Fig. [TOl 

A stringent test of the computational framework to de- 
scribe dipolar molecules is provided by the analytic result 
/icr = 0.639 eao for the fixed dipole (/ — oo) [40|. To this 
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FIG. 10. Similar to Fig.[8]but for the second excited J" = O3 
state of LiCl^ . 



end, we performed BEM calculations for a dipolar system 
at steadily decreasing moments of inertia [H, [l^. For 
each value of /, the dipolar anion energies have been cal- 
culated for 1080 values of p. in the interval 0.6 < /i < 3.0. 
Only 170 energies satisfying the subthreshold condition 
E < Eiim = — 10~®Ry were retained to minimize the 
numerical error. These energies correspond to an inter- 
val A/i ~ 0.377 of the dipole moment. We checked, that 
in this energy interval, /icr can be obtained by using the 
expression 



(41) 



to extrapolate the calculated energy down to = 0. One 
should stress however, that an excellent energy fit in the 
subthreshold region does not guarantee an excellent es- 
timate of the critical dipole moment. The values of ficr 
extracted by this extrapolation procedure can be consid- 
ered reliable only if A/i, which depends on the chosen 
precision Eum, is close to the critical dipole moment. In 
the cases studied, this criterion is approximately satisfied 
only for the ground state and the first excited O'^ state. 
The critical dipole moments for these states in anions 
with the dipole length s = 4ao are shown in Table U for 
various moments of inertia. The agreement with the ana- 
lytic limit is excellent for the ground state configuration, 
and is fairly good for the first excited 0+ state. This is 
very encouraging, considering the slow convergence with 
/ and various sources of numerical errors in the E 
regime p^ . 



VIII. RESULTS FOR SPECTRA AND RADII OF 
DIPOLAR ANIONS 

Energies and r.m.s. radii of the lowest bound 0+ and 
1~ states of Lil~, LiCl~, LiF~, and LiH^ dipolar anions 
predicted in this study are listed in Table [ITl 



TABLE I. Critical dipole moments for dipolar anions in the 
two lowest 0"*" states calculated in this work (BEM) and in 
Ref. [3 for the charge separation s = 4ao and different mo- 
ments of inertia I. The analytic result at J — > co j^, is 
ficr = 0.639 eao. 



I{meal) 


BEM 


(eao) 
Ref. [19] 


per 

BEM 


(eao) 
Ref. [19] 


10^ 


0.937 


0.843 


1.024 


1.515 


lO'' 


0.674 


0.750 


0.633 


1.145 


10« 


0.639 


0.715 


0.622 


0.974 


10^° 


0.639 




0.622 




10^5 


0.639 




0.62 





TABLE 11. Energies and r.m.s. radii for 0^ and 1~ bound 
states of selected dipolar anions obtained in DIM {Imax ~ 9) 
and BEM. The parameters of the calculation are given in Sec. 
IVII The numbers in parentheses denote powers of 10. 



Anion 


state 


E{Ry) 
DIM BEM 


DIM 


(ao) 
BEM 


Lir 


0^ 


-5.079( 


-2) 


-5.023( 


-2) 


7.569(0) 


7.620(0) 






-9.374( 


-4) 


-1.037( 


-3) 


5.112(1) 


4.759(1) 




ot 


-1.502( 


-5) 


-1.797( 


-5) 


3.719(2) 


3.308(2) 






-5.079( 


-2) 


-4.995( 


-2) 


7.569(0) 


7.641(0) 






-9.291( 


-4) 


-1.023( 


-3) 


5.112(1) 


4.886(1) 




I3" 


-1.261( 


-7) 


-1.099( 


-5) 


3.423(3) 


3.464(2) 


Licr 




-4.483( 


-2) 


-4.483( 


-2) 


7.885(0) 


7.894(0) 






-7.374( 


-4) 


-8.241( 


-4) 


5.632(1) 


5.017(1) 




ot 


-7.051( 


-6) 


-9.907( 


-6) 


5.124(2) 


4.106(2) 






-4.482( 


-2) 


-4.458( 


-2) 


7.885(0) 


7.915(0) 






-7.241( 


-4) 


-8.067( 


-4) 


5.633(1) 


5.337(1) 






-3.062( 


-7) 


-8.159( 


-7) 


2.066(3) 


8.831(2) 


LiF" 


0^ 


-2.795( 


-2) 


-2.983( 


-2) 


9.117(0) 


8.991(0) 




0^ 


-3.022( 


-4) 


-3.525( 


-4) 


8.098(1) 


7.501(1) 




0^ 






-6.101( 


-8) 




3.363(3) 




ir 


-2.793( 


-2) 


-2.968( 


-2) 


9.117(0) 


9.010(0) 




I2" 


-2.782( 


-4) 


-3.277( 


-4) 


8.124(1) 


7.520(1) 


LiH" 




-2.149( 


-2) 


-2.370( 


-2) 


1.011(1) 


9.698(0) 






-1.491( 


-4) 


-1.922( 


-4) 


1.058(2) 


9.297(1) 




ir 


-2.142( 


-2) 


-2.353( 


-2) 


1.011(1) 


9.717(0) 




^2 


-7.942( 


-5) 


-1.231( 


-4) 


1.146(2) 


9.591(1) 



One can see that for each total angular momentum 
there are at most three bound eigenstates in each 
system. The r.m.s. radius of an electron cloud shows a 
spectacular increase with decreasing the binding energy 
of the state. For the subthreshold states, such as 0;^ and 
I3", the radius is of the order of hundreds to thousands 

Energy spectra and radii of dipolar anions do not 
change significantly in the limit imax ^ 00. Usually, 
the extrapolated results for both E and r^ms agree very 
well with those in Table |lT] {imax = 9). For instance, 
the extrapolated values for the state in LiH~ are 
E = -7.931 • 10-5 Ry and r.^s = 1-147 • 10^ oq. 
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The DIM and BEM results are generally consistent for 
both energy and radii though significant quantitative dif- 
ferences persist for excited, weakly-bound states of anions 
where the DIM is not expected to work. In the case of 
LiF~, the BEM predicts the existence of the third 0;^ 
state at an energy —6.1 • 10^* Ry, which is absent in 
DIM. 

It is instructive to compare our DIM results with those 
found in Ref. [l^] using a similar approach. Tablc HlIl lists 
energies of the lowest 0+ bound states of Lil~, LiCl", 
LiF~, and LiH~ dipolar anions obtained in both studies, 
and Table |IV] shows the adopted values of dipole mo- 
ments. 



TABLE III. Energies for 0^ bound states of selected dipo- 
lar anions obtained in DIM in this work {imax = 9) and in 
Ref. (23 |. The numbers in parentheses denote powers of 10. 



Anion 


state 


E (Ry) 


Ref. [22] 


This work 


Lil" 




-5.079(-2) 


-4.998(-2) 




ot 


-9.374(-4) 


-1.022(-3) 




ot 


-1.502(-5) 


-1.999(-5) 


Licr 


ot 


-4.483(-2) 


-4.483(-2) 




ot 


-7.374(-4) 


-7.497(-4) 




ot 


-7.051(-6) 


-9.775(-6) 


LiP- 


ot 


-2.795(-2) 


-2.793(-2) 




ot 


-3.022(-4) 


-3.366(-4) 




ot 




-8.746(-7) 


LiH" 


ot 


-2.149(-2) 


-2.352(-2) 




ot 


-1.491(-4) 


-1.926(-4) 



The two calculations agree reasonably well for the 
lowest-lying states; some difference stems from slightly 
different dipole moments used in Ref. [l^ and here. In- 
deed, while the charge separation in both studies was ad- 
justed to reproduce the experimental ground state energy 
of LiCl~ , the fitted values of s in both calculations are 
different: s = 0.3335 ag in Ref. ^2^] and sdim = 0.336 qq 
here. 

The largest deviations, seen for weakly-bound states, 
can be traced back to the cutoff value of the electron 
orbital angular momentum when solving CC equations. 
In Ref. [22], adopted imax was small, typically Imax — 4 
[3^, whereas it is fairly large, imax — 9, in our work. 
As seen in Fig. [S] and discussed in Sec. IVII[ energies of 
weakly-bound states obtained in DIM do converge slowly 
with imax- Therefore, calculations employing low imax 
values cannot be useful when performing extrapolation 



imax ^ 

IX. CONCLUSIONS 

In this study, we applied the theoretical open-system 
framework based on the Berggren ensemble to a problem 
of weakly-bound states of dipole-bound negative anions. 



TABLE IV. Dipole moments of selected dipolar anions 
adopted in this work and in Ref. [23 |. 



Anion 


This work 


fi (eao) 


Ref. [22] 


Lir 


2.911384272 




2.911384272 


Licr 


2.805158089 




2.793355179 


LiP- 


2.472316049 




2.478610934 


LiH^ 


2.313370205 




2.321238811 



The method has been benchmarked by using the tradi- 
tional technique of direct integration of CC equations. 
While a fairly good agreement between the two methods 
has been found for well-bound states, the direct inte- 
gration technique breaks down for weakly-bound states 
with energies \E\ < 10~*Ry, which is comparable with 
the rotational energy of the anion. For those subthresh- 
old configurations, the Berggren expansion is an obvious 
tool of choice. 

The inherent problem of the DIM is the lack of stabil- 
ity of results when the number of channels is increased. 
Indeed, the method breaks down when the channel or- 
bital angular momentum becomes large; around ic — 10. 
This can be traced back to the applied matching con- 
dition. We demonstrated that this pathology is absent 
in BEM. Here, the rapid converge with i is guaranteed 
by an effective softening of the interaction through the 
momentum cutoff kmax, which suppresses contributions 
from high-^ partial waves. 

The future applications of BET will include the struc- 
ture of quadrupole-bound anions [41I, [90l - l92j and the con- 
tinuum structure of anions, including the characteriza- 
tion of low-lying resonances. 
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